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Abstract 

The bypass duct of an aircraft engine is a low-pass 
filter allowing some spinning modes to radiate out- 
side the duct. The knowledge of the radiated modes 
can help in noise reduction, as well as the diagno- 
sis of noise generation mechanisms inside the duct. 
We propose a nonintrusive technique using a circular 
microphone array outside the engine measuring the 
complex noise spectrum on an arc of a circle. The ar- 
ray is placed at various axial distances from the inlet 
or the exhaust of the engine. L’sing a model of noise 
radiation from the duct, an overdetermined system 
of linear equations is constructed for the complex 
amplitudes of the radial modes for a fixed circum- 
ferential mode. This system of linear equations is 
generally singular, indicating that the problem is ill- 
posed. Tikhonov regularization is employed to solve 
this system of equations for the unknown ampli- 
tudes of the radiated modes. An application of our 
mode detection technique using measured acoustic 
data from a circular microphone array is presented. 
We show that this technique can reliably detect ra- 
diated modes with the possible exception of modes 
very close to cut-off. 
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List of Symbols 


a 

radius of circular microphone array 

A(w,0') 

a complex function defined in equa- 
tion (17) 

-1, 

the vector (A q i , A ql , . . . , A qn ) T 

4 

amplitude of mode (m, n), rn — aB t 

qV 

B 

fan blade number 

BPF 

blade passage frequency (BU/2tt) 

c 

speed of sound 

C(m t n, c) 

t a parameter defined in equation (12) 

(iq 

the vector 


a paraniter defined in equation (16) 

f 

frequency 

G 

a j x n matrix, j > n, see equation 
(23) 

G* 

the Hermit ian adjoint of G 

G 

the linear operator mapping X into 
Y whose matrix is G 

Ovw 

element of the matrix G 


Bessel function of first kind and or- 
der rn 

k 

wavenumb er , a BQ/c 

k a (m,n) 

axial wavenumber of mode (rn, n), 
see equation (22) 

k r (m,n) 

radial wavenumber of mode (m,n), 
see equation (19) 

m 

circumferential mode number 

M 

duct Mach number, M < 0 for inlet 
radiation 

n 

radial mode number, n = 1,2,... 

Ng 

the nullspace of linear operator G 

p'(x, t) 

acoustic pressure 


temporal Fourier transform of 
p'(x,t) 

P{x,aBQ] 

|s ee equation (14) 

A 

Fourier transform of P in azimuthal 
variable 0 r , see equation (15) 

q 

a parameter defining circumferential 
mode m in the relation m = aB + 
q\\ g — 0, ±1, ±2, . . , 
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r 

radial coordinate of a point on the 
duct inlet or exhaust surface, see fig- 
ure (2) 

ro 

duct radius 

Rg 

the range of linear operator G 

Rq 

distance of a microphone from the 
duct inlet or exhaust center 

R 

distance between a microphone and 
a point on the the inlet or exhaust 
surface 

ret 

as a subscript, retarded time 

s 

the dimension of the range of G 

S 

in dS , surface area 

SYD 

Singular Value Decomposition 

t 

time 

T 

as a subscript, Tikhonov 

u 

axial velocity of the sound wave 

Ui 

basis vectors in the sub space X \ 
V g ) i — 1 , 2 , . . . , 

Vi 

basis vectors in Re , i = 1, 2, . . . , .5 

V 

number of vanes in stator, number 
of circumferentially symmetric dis- 
turbances interacting with the rotor 
blades 

X 

axial distance from inlet or exhaust 
center 

X 

the ri- dimensional linear space where 
the solution vector A q is sought 

Y 

the j-dimensional linear space where 
the vector d q lies 

Greek Symbols: 

a 

multiple of BPF 

C\T 

Tikhonov parameter 

0 

a/1 - A/ 2 

9 

angular coordinate of a point on the 
inlet or exhaust surface in cylindrical 
polar coordinates, see figure (2) 

9' 

angular coordinate of a microphone 
in cylindrical polar coordinates, see 
figure (2) 

£(m,n) 

cut-off ratio, sec equation (20) 

/‘i 

singular values of matrix G, — 

1 , 2 ,..., n, pf arc the eigenvalues of 
G*G, & > 0 

A) 

density of air in the exterior region 
of the duct 

0 

angle that the line joining a micro- 
phone to the center of the duct inlet 
or exhaust makes with the or-axis 

n 

angular velocity of the rotor 

Other Symbols: 

INI 

Li , L‘ 2 > or norm of a vector 


■ \ • as in X \ Nn , set subtraction 

< . > complex inner product. 

1. Introduction 

The detection of the radiated modes from an air- 
craft engine can help in the diagnosis of noise gener- 
ation mechanisms inside the duct and provide guid- 
ance for noise control. A duct is a low-pass filter 
that allows some of the spinning modes inside the 
duct to escape outside. In general, each radiated 
mode has a specific directivity pattern. The com- 
bined interference of all the radiated modes results 
in a complicated spatial pattern that rotates at the 
engine speed. The problem of the discovery of any 
information about the radiated modes is known as 
an inverse problem. The study of inverse problems is 
currently an active area of research. These problems 
appear in many areas of science, engineering, and 
medicine. For example, magnetic resonance imag- 
ing (MRI) in medicine, which allows physicians to 
obtain images of the internal organs of patients, can 
be classified as an inverse problem. Other examples 
appear in oil and mineral prospecting, seismology, 
and satellite imaging. 

One method of mode detection was proposed by 
Pickett et al. in the mid-Seventies. 1 It uses a micro- 
phone rake rotating at steady speed in the duct in- 
let or exhaust. This method was successfully imple- 
mented by Heidelberg and coworkers who designed 
a microphone rake rotating at a fraction of the en- 
gine speed. 2 The rotating rig size is dependent on 
engine size, requiring sophisticated mechanical and 
electronic components. In addition, the wake from 
the rake may interact with fan blades if used in the 
inlet. The resulting interaction noise can be sepa- 
rated from the engine noise in most cases, but may 
complicate the interpretation of the results. We 
present another method which employs a circular 
microphone array outside the engine and. thus, is 
nonintrusive. The associated measurement rig can 
be used for a range of engine sizes. In addition, the 
microphone rig is very simple, using some stepping 
motors to move the microphones (located on a cir- 
cular arc) circumferentially and axially. 

We employ Rayleigh's piston in a wall formula to 
model the radiation from the inlet or exhaust plane 
of the duct. For each multiple of the blade pas- 
sage frequency (BPF), the complex temporal Fourier 
transform of the acoustic pressure is first obtained 
over an arc of a circle at a fixed axial distance from 
the duct lip. The Fourier decomposition in the az- 
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imuthal direction of the temporal transform can then 
be related to the amplitudes of the radiated modes 
at the duct inlet or exhaust. For a given circumfer- 
ential m, a single linear equation is obtained in terms 
of the amplitudes of the radial modes which propa- 
gate outside the duct. Repeating this calculation for 
several axial distances, a set of simultaneous linear 
equations for the amplitudes of the radial modes for 
the circumferential mode rn is obtained. In princi- 
ple, this system of linear equations can be solved to 
get the amplitudes of the radial modes. There are, 
however, some problems. For a given circumferential 
mode, the number of radial modes that will propa- 
gate outside of the duct is not known. Therefore, 
we propose to take more axial measurements than 
is necessary (i.e. we work with an overdetermined 
system of equations). Furthermore, the problem is 
ill-posed and techniques for extracting the proper 
information from the linear equations must be em- 
ployed. We use the Tikhonov regularization tech- 
nique which is a well known method of dealing with 
ill-posed problems. 3, 1 

In Section 2, we present the detailed mathemati- 
cal derivation behind the mode detection technique. 
There are practical problems in finding the phase 
of the complex temporal Fourier transform of the 
acoustic field for a multiple of the BPF because we 
can not measure sound simultaneously at ail points 
of the space. This necessitates the use of a reference 
microphone and the cross spectrum to extract the 
phase information. In Section 3, we discuss the ill- 
posedness of the problem and how we remedy this 
problem by regularization. Results from an experi- 
ment on a 12 inch diameter model are presented in 
Section 4. The calculations show that our mode de- 
tection method works well for modes that are not 
very close to the cut-off condition. Finally, conclud- 
ing remarks follow in Section 5. 


The number of blades of the rotor is B and H is 
the shaft angular frequency. The blade passage fre- 
quency (BPF) is. therefore, BQ/2n. The number 
of stator vanes is V', which can also be assumed to 
be the number of axially symmetric disturbances in- 
teracting with the rotor. The circumferential mode 
numbers, m. due to rotor-stator interaction are given 
by the relation 

m^aB + qV q~ 0, ±1,±2,... (1) 

where a is the multiple of BPF. For a given a, we 
note that q specifies m, i.e. m — m(q). In what 
follows, we find that it is convenient to work with 
q rather than m. Our aim is to obtain a system of 
simultaneous linear equations for the radial modes 
associated with a fixed o and q that radiate outside' 
the duct. Although we refer to duct inlet noise ra- 
diation, the analysis is also valid for duct exhaust 
with the appropriate duct Mach number, A/, axial 
wave number, k a , and radial wave number, k r . 

We use Rayleigh's piston in the wall formula as 
the model of noise radiation from the inlet 

4 V(M) = [ dS (2) 

'/inlet area + *■ 

where u(r, 0,0, £) is the fluctuating axial velocity of 
the sound wave at the inlet, x — 0, and R is the 
radiation distance as shown in figure 2. The origin 
of the cylindrical polar coordinate (r,0,.r) is at the 
center of the duct inlet with the microphone coordi- 
nates at (a ) 6 ( ,x). The acoustic pressure in the duct 
is described by 

= e ickB ( Qt -°)x 

(3) 

E, - Vm [Mm, n)r] e -.-[*vw.(m,»)xl 


2, Mode Detection Theory 

In this section, the mathematical theory of the 
mode detection, which was originally discussed by 
Farassat and Myers 5 is presented. Figure 1 shows 
the circular microphone array (CMA) of radius a at 
a distance x from the duct, radius ro, inlet (or ex- 
haust). The microphone array plane is parallel to 
the duct inlet plane and the center of the array cir- 
cle lies on the extension of the duct axis as shown in 
figure 1. It is assumed that we have a uniform flow 
with the Mach number M inside the duct. Since the 
positive .r-axis points out of the duct, M < 0 for 
the inlet and M > 0 for the exhaust noise study. 


where J m (-) is the Bessel function of the first kind 
of order m, k a is the axial wave number, and k r is 
the radial wave number. Note that both k a and k r 
are functions of the circumferential mode number m 
and radial mode number n. In equation (3), A qn is 
the amplitude of the mode (m,n) ( i.e. A qn stands 
for A mri ). We are summing over q in this equation 
and m ~ aB + qV = r n(q) everywhere. 


The relation between j/(m,n) and u(m,n) is 
found from the linearized momentum equation as 


u(m,n) 


k a (m, n)p'(m, n) 


( 4 ) 
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Using this equation and equation (3), the axial ve- 
locity of sound at the inlet, x — 0, is 

T, n Y,q A qnka{m,n)J m [k a (m,n)r]e- iqV9 
which simply gives 

tx(r, 9, 0, £) — i<xBSlu{r } 9 } 0, t) (6) 

Next, we assume that the duct radius ro satisfies 
r 0 << Ro where i?o is the distance from a micro- 
phone on the array to the duct inlet center. Then, 
since r < r 0 , we can show that 

R « Rq - r sin ip cos (9 -- 9 f ) (7) 

where sin */? = a/Ro and ip is the angle that the line 
joining a microphone to the duct inlet center makes 
with the ;r- axis (see figure 2). From equations (5-7), 
we obtain 


where 

. f r ° _ .a Bttr . ri / m , 

C(m, n , ip) = / rJ m (— sin ip)J m dr 

Jo C 

( 12 ) 

Our mode detection technique is based on equa- 
tion (11) which we will discuss below. 

The complex temporal Fourier transform of 
p r (x, t ) at a frequency of a times BPF is 

E„ E, C(m,n, mMm, „) e -*W -»/*> 

(13) 

Let us now define 
P(x, tvRfi) 

= -2iRoe' aB(nRo/cf0 ''’ r/2) p'(x,QBn) (14) 


[u] = 


Aj 


n q 

i[fcr sin tp cos (9-0')- m(9—0')] 


r) x 


(8) 


where k = aBQ/c and fc a and are functions of 
(m, n). 

We can write the surface integral in equation (2) 
as follows: 


/ 

j in 


. . .dS = 

inlet area JQ 


prq r2n 

/ rdr 
Jo Jo 


de (9) 


Using the approximation R « Rq in the denominator 
of equation (2), the only term in [*i] ret dependent 
on the angle 9 is the last exponential function in 
the summation in equation (8). This term can be 
integrated analytically in the variable 9 as follows: 


r 2K 

Jo ■ 


sin V’ cos(0- 0 f ) - m(0 — 


2i re imn /~J m (kr sin ip) 


(10) 


The acoustic pressure at the microphone location 
x = (a } 9 f ,x) is thus given, from equation (2), as 


~Ro/c)-0' + ir/2] y 

E n E,C’(tn I n,^).4, n * 0 (m,fi)e- < * v l*'-'/ 3 J 


(ID 


If we find P(x,aBQ) for 0 < 9' < tt, then the 
complex Fourier Transform of P in the azimuthal 
direction, 0 f gives 

A = - r P(x,aB{l)e iqre ' d$' (15) 

^ Jo 

- e' qVw/2 [C(m, n, rf>)k a (m , n)] .4„„. 

n 

From this result, we find 

A, = e- iqVK/2 P v (16) 

= P(m, n, ip)k u (m, n)] .4,,,. 

n 

This is a linear equation for the amplitudes of all 
radial modes associated with a fixed circumferential 
mode, 7?7 = aB + q\\ that radiate outside the duct. 
Each axial position of the array supplies one linear 
equation for A nq . Measuring the noise radiated from 
the duct at many axial locations gives an overdetcr- 
mined system of linear equations for the known 
( q fixed) for the associated radial modes. One then 
solves this system of linear equations to get the ra- 
dial mode amplitudes A qn . 


Remarks 


1. Equation (11) shows that the acoustic field of 
the noise radiated from the duct may be represented 
in the far field as 

n - ik Rq 

p'(x,t) = Mip,0')e iaB{Ut ° } (17) 
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where A {ip,0') is a complex function which is pe- 
riodic in the azimuthal angle, 9 r with periodicity 
2n/Y. This result indicates that the information 
about A qn is hidden in A (0, 9'). Our mode detection 
technique depends on obtaining one linear equation 
for A qn (for fixed values of q and ip) from the Fourier 
transform of A (ip ) 9 t ) in the azimuthal direction 9 r . 


3. Ili-Posedness and Remedy 

We are finding the complex amplitudes of the ra- 
diated modes A qn by solving an overdetermined sys- 
tem of equations 

GA q = d q (23) 


2. Note that the coefficients of the unknown am- 
plitudes A qn in equation (16), C(m, n, ip)k a {m } n), 
are known in advance and numerically given by q. n, 
and v. These coefficients do not depend on experi- 
mental data. The measured data are used in finding 
D q which is obtained schematically as follows 

— complex const, x e ~ tq}> K ^ 2 x 


where G is a j x n matrix, j > n, 
Aq — (A^i , Aq 2 ? ■ ■ * j A^n) and d q 
(Dq(tp\) } • • ■ , Dq(ipj)). Equation (23) is the 
matrix form of the system of linear equations based 
on equation (16) and evaluated at axial distances 
j*i , ^ 2 , . . . , xj corresponding to ip — ^’i , 02 1 ■ • • > 4’j 
such that tpi > 02 > - -- > 4 J j ■ We note that if the 
elements of the matrix G are denoted by g vw , then 


{spatial F.T. [temporal F.T. of p'(x, £)]} 

(18) 

3. In general, one must have some knowledge of 

V in advance. For example, if wc have stator vanes, 
then V can be taken as the number of stator vanes. 
Otherwise, one can only guess about circumferen- 
tially symmetric disturbances that generate a de- 
tected mode. This is not necessarily a weakness of 
our method because the detected modes can guide 
the user to search for possible circumferential dis- 
turbances interacting with rotor blades both inside 
and outside the duct. For example, the ribs on an 
inflow control device intended to break down the at- 
mospheric turbulence may produce wakes that can 
interact with rotor blades in the duct. In this case, if 
the ribs are placed at equal angles circumferentially, 

V should be taken as the number of ribs. 

4. The radial wave number k r (m,n) is found from 
the solution of 


J’ m [A- r (m,n)r 0 ] - 0. (19) 


We index n starting with n = 1 rather than n = 0. 
The cut-off ratio, £(m, n) of the mode (m, n) is given 


by 


£(m, n) 


k 

pk r (m } n) 


( 20 ) 


where p = V 1 - A/ 2 . The axial wave number. 
k a (m,n ) , is found from 


M m ,n) - 


~M ± y/\ l/£ 2 (m,n)| . (21) 


For inlet radiation and the .r-axis pointing outward, 
since M < 0 and k a {m,n) > 0 the relation becomes 


*.(m, n) = A [| M | +y/\ - l/e(m,n)] . (22) 


Qvw — (m , IV , ip v ) k a (tti , IV ) (2d) 


where C(m } w,ipv) and k a (m t w) are given by equa- 
tions (12) and (22) with appropriate change of no- 
tation. We will show below that g lw can be written 
in a particularly simple form which can be useful 
in applications. Using equation (5.54), page 634 of 
Gradshteyn and Ryzhik, 6 we have 

C(m % w } y> v ) - 9 1 [3 W r 0 J m {a v r 0 ) J m _ i (p u ro ) 

ro«/m — i (f*vf*o)«An (/^u'^o)] (25) 

where 

a v — sin ip v = k sin xp v (26a) 

c 

Pw = k r (m,w). (26b) 


Now we use the following relation (eq. 8.472, p.967, 
Gradshteyn and Ryzhik 0 ) 

zJ,n{z) = zj-j m {z) + mJ m (z) (27) 
(t~ 


and the fact that J r m (p w r o) = 0, because of the 
boundary condition, to simplify the right side of 
equation (25) and obtain 


C{m,nyi\P) 


Jm{Pw1* o) r J ( \ 

2 q2 m (^V^o) 

a v ~ Pw 

otcro ./ m (a ( ,r 0 )] 

J m V J" o)J m ( /3 it >'o ) 


(28) 


The expression on the right is not singular when 
a v -> p a because then J r m (P w ^o) ~ 0- We, there- 
fore, have 


gvw = 


K ( m , n) J r m (a t; r 0 ) J m {p w r 0 ) 
3 2 a 2 


(29) 


O 


( 31 ) 


giving the simplest form in which g vw may be writ- 
ten. 

In general, we do not know whether a radial mode 
that can propagate in the duct was actually gener- 
ated. Therefore, we must have the corresponding 
column in matrix G. We propose taking more mea- 
surements in the axial direction than is necessary to 
find all of the radial modes. Because of measurement 
and numerical errors, the overdetermined system of 
equations (23) is not consistent. It is a good idea 
to keep all the linear equations and use a method to 
extract the amplitudes of the radiated modes from 
the overdetermined system of equations (23). Even 
with very careful measurements, it is not easy to get 
precise information on mode amplitude because the 
matrix G has zero or near zero singular values and, 
therefore, is ill-conditioned. 

To obtain a meaningful solution for A q from 
equation (23), we must use a regularization tech- 
nique such as singular value decomposition (S\ D) or 
Tikhonov regularization. 3,1 Both methods involve 
first finding the n singular values fx\ > fx\ > . . . > 
fi 2 n > 0 which arc the eigenvalues of the n x n matrix 
G*G. Here, G* is the Hermitian adjoint of G. It is 
possible that the singular values beyond some index 
are all zero or very small. When the singular val- 
ues are very small, any errors in the elements of d q 
are magnified greatly, resulting in large errors in the 
solution of equation (23). This situation is very com- 
mon in practice, a point emphasized by Tikhonov in 
numerous publications. 7 Tikhonov suggested a reg- 
ularization technique known by his name which we 
have used in our work. 

Referring to figure 3, let A" be an n-dimensional 
space and Y a j -dimensional space. The matrix G 
represents a linear operator G : X — > Y. In general, 
the nullspace of G, No, is not empty and the range, 
Rq , of G is not the entire space I’, *.e., Re / Y. 
It is entirely possible that d q in equation (23) is not 
in but has a projection in Rq- We know that 
X \ Nq — > Rq is an isomorphism and, therefore, 
X\Nq and Re have the same dimension s < n. The 
singular values //[ . //•_>, .... j.i s are all non-zero. Both 
SYD and Tikhonov regularization depend on the 
construction of an orthonormal basis (rq , U 2 , . ■ • , u s ) 
in X \ Nq and another (ci , t ?2 , . - - , v s ) in /?o* Then 
in the method of SYD, the solution of equation (23) 
is interpreted as the solution of 


G*G.4, - G*d q (30) 


which always exists. It is given by 
5 ^ 1 

(^).svd = 7T ^ v i > Uj 

i=l to 

where < . > is the inner product. 


It is possible that some of the singular values be- 
yond some index are very small and then small errors 
in d q are magnified resulting in large errors in A fJ . 
Tikhonov proposed to use the following result as the 
solution of equation (23) 


(^<i) T ~ Y. 


i= 1 


Vi 

otr + vi 


< d q , Vj. > u t (32) 


where ar > 0 is a constant known as the Tikhonov 
parameter. This parameter is usually of the order of 
the square of the smallest nonzero singular values. 
It can be shown 3 that (A q )r minimizes the function 


||G(.4 ? )t - d q \\y + ot||(*4Jt||.y (33) 

which means that we are controlling the error 
|| G(A q )r — d q \\ y and the norm of the solution 
IK*4 (/ )t||.y* There are many articles discussing the 
applications of the Tikhonov regularization method 
and the process of selecting the parameter ar. 4 We 
have chosen this method for the mode detection 
technique. 


4. Some Results 

We now present some results from an experiment 
on a 16-bladed fan with 20 vanes in a duct of radius 
r 0 = 0.1524 m (6 inches). The inlet flow Mach num- 
ber was 0.605. The rotor speed was 16,900 RPM 
corresponding to BPF = 4507 hz. The radius of the 
circular microphone array was 0.912 m (3 ft.) with 
18 microphones located at 10 degree intervals. The 
microphones rotated in steps of 0.5 degrees covering 
10 degrees for a total of 360 points on the arc of half 
a circle. Ten axial locations wore used in this test, 
the set-up of which is shown in figure 4. 

In order to measure the phase of the complex tem- 
poral Fourier transform, a stationary reference mi- 
crophone is needed to make simultaneous noise mea- 
surement with microphones on the array. The cross- 
spectral analysis of two microphones (reference and 
an array microphone) were used to get the phase in- 
formation of the temporal Fourier transform. It is 
not at all obvious where to place the reference mi- 
crophone. One would like to place the reference mi- 
crophone in the peak directivity position of a mode. 


6 



This is not, in general, possible because we are look- 
ing for many modes with different directivity pat- 
terns. For this test, we have used three reference 
microphones at xp = 30, 50, 70 degrees, as shown 
in figure 4, at a constant distance of 1.83 m (6 ft) 
from the fan axis. Theoretically, the three reference 
microphones should give the same mode amplitudes 
for any radiated mode. 

We were interested in the rotor-stator interaction 
modes for BPF, 2BPF, and 3BPF. The RPM of the 
fan rotor was varying within 0.6 percent of the spec- 
ified value. For modes that radiate outside the duct, 
( i.e . the modes that are not very close to cut-off), 
we should be able to find their amplitudes with this 
small variation in RPM. The reason is apparent from 
equation (17). The temporal Fourier transform of 
gives a complex quantity independent of ft. 
Thus, small variations in ft should not introduce er- 
rors in the amplitudes of the radiated modes. 

Test Results 

1. The phase of the temporal complex Fourier 
transform is always given relative to the reference 
microphone. Therefore, for a fixed axial location 
of the array and over the 360 microphone positions 
on the half-circle, the phase distributions from two 
reference microphones should show a constant shift 
between them. This is a good test of the phase mea- 
surement technique and instrumentation. For BPF 
and 2BPF, the shift varied by approximately 5 de- 
grees. For 3BPF, the shift variation was larger at 
approximately 10 degrees. This order of accuracy 
does not seem to be satisfactory. This can be a 
source of error in finding the mode amplitudes and 
the disagreement of the results from different refer- 
ence microphones. 

2. The following are the possible interaction 
modes that can radiate noise outside the inlet: 

BPF: (-4,1) to (-4,3) 

2BPF: (-28,1), (-8,1) to (-8,7), (12,1) to (12,5) 

3BPF: (-32,1) to (-32,3), (-12,1) to (-12,10), (8,1) 
to (8.12), (28,1) to (28,4) 

We will present results for some of the modes shown 
here. Since we had 10 axial stations in the test, we 
cannot compute the amplitudes of modes (8,11) and 
(8.12) for 3BPF. 


3. Ill-posedness of the problem becomes obvious 
if we look at the singular values of the matrix G 
(square roots of the eigenvalues of G*G). For cir- 
cumferential mode -12 at 3BPF, these are: 


n 

Pi 

1 

4.42o-4 

2 

3.63c-4 

3 

2.98o-4 

4 

2.37o-4 

5 

2.10e-4 

6 

1.73e-4 

7 

5.58e-5 

8 

9.89o-7 

9 

3.440-9 

10 

1.52e-12 


Table 1. Singular values for radial modes of m = -12 
at 3BPF 


It is quite obvious that we have a very ill-conditioned 
matrix G*G. Similar behavior has been observed 
for other modes. This is a clear indication of the 
ill-posedness of the mode detection problem. 

The general rule for using a regularization tech- 
nique is that one must have a good idea about the 
errors in equation (23). In our case, we must be 
concerned with errors in the measured quantity, d q < 
because the elements of G are known analytically 
to machine accuracy. A phase error of A <f> in the 
temporal Fourier transform introduces an error in 
d q that is magnified when finding the mode ampli- 
tude by a factor 1//0: if the error vector lies along 
the basis vector, Vi, in Y space (see equation (31)). 
Because the singular value can be very small, the 
error introduced in the computed mode amplitude 
can be relatively large. 

After testing various values of the Tikhonov pa- 
rameter, otj it was decided that a value of ot — 
10” 10 was a good choice, essentially affecting /o ~ 
10 ” 5 and smaller. The following results have been 
calculated using this value of the Tikhonov parame- 
ter. 

4. Mode amplitudes for BPF from the three ref- 
erence microphones are shown in figure 5. In this 
figure, and those following, we will not show the am- 
plitude of the radial mode nearest to cut-off (highest 
radial mode) unless there is only one mode that is 
cut on. Figure 5 shows the good agreement between 
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the three reference microphones in determining the 
mode amplitudes. 

Figure 6 shows amplitudes of the radiated modes 
for 2BPF and figure 7 shows similar results for 
3BPF. In both of these figures, we see that there 
are some large amplitudes for modes near cut-off. 
Because these modes have predominantly in-plane 
directivity, these large amplitudes are almost defi- 
nitely caused by the ill-conditioning problem of the 
matrix G. Otherwise, there is general agreement be- 
tween the amplitudes of the modes computed from 
the three reference microphones, although there arc 
some notable disagreements. It appears that errors 
in the phase measurements of the temporal Fourier 
transform of p'(x, t) are the main cause of the dis- 
crepancies. 

On the whole, we can claim that the mode detec- 
tion technique discussed in this paper appears to be 
doing a good job of giving us the mode amplitudes. 
Further refinements of the technique are possible as 
discussed below. 

5. Concluding Remarks 

We have presented a mode detection theory for 
noise radiated from the inlet or the exhaust of an 
engine. We have shown that the spatial Fourier 
transform in the azimuthal direction of the complex 
temporal Fourier transform of the acoustic pressure 
p'(x, t) has the information about the amplitudes of 
all the modes radiated outside of the duct. The 
problem of mode detection is ill-posed and a reg- 
ularization technique such as SVD or the Tikhonov 
method must be used to obtain meaningful results. 
In using any regularization technique, one must have 
an estimate of the errors in the measured quantities 
used as input to the linear system whose solution 
gives the amplitudes of the radiated modes. This 
estimate is not always available resulting in some 
loss of accuracy. 

The example given in this paper shows that our 
mode detection technique is capable of finding the 
amplitudes of the radiated modes. It can be a very 
useful tool in the control and reduction of engine 
noise. There are many areas of improvement to make 
this method more useful such as: 


i. Improving the accuracy of phase measurements 
of the complex temporal Fourier transform of 
the noise at microphone locations. 


ii. Changing the duct radiation model to include 
the effect of refraction of sound by the inflow 
or exhaust stream. This will complicate the al- 
gebra for setting up the linear system of equa- 
tion (23), but will undoubtedly be a more real- 
istic model than the one currently employed. 

iii. Exploring other regularization techniques. 3,1 
Currently, many problems of science and tech- 
nology are ill-posed and many publications ad- 
dress the selection of the regularization method. 
Some of these are good candidates for the mode 
detection problem. 
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Figure 1 . Microphone array and duct geometry 
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Figure 2. Description of symbols used for mode detection theory 
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r * > Isomorphism 

y\n g = r g -^Vx\n g = r g *J 

DIM R g = DIM R G * = s < n 


Figure 3. The behavior of mapping G and G* between linear spaces X and V , Nq 
of G, Nq is nullspace of G * , Rg is range of G 
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Figure 4. Schematic of the mode detection rig 
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Figure 3. Radial Mode Amplitudes for BPF using the three reference microphones 
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Figure 5. Radial Mode Amplitudes for BPF using the three reference microphones (cont.) 
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Figure 6. Radial Mode Amplitudes for 2BPF using the three reference microphones (cont.) 
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Figure 7. Radial Mode Amplitudes for 3BPF using the three reference microphones (cont.) 
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